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Abstract 

Starting from fixed-order perturbation theory (FOPT) we derive expres- 
sions for the heavy-flavour components of the deep-inelastic structure func- 
tions (Fi } H{x,Q 2 ,m 2 H ),i = 2,L;H = c,b,t) in the variable-flavour number 
scheme (VFNS). These expressions are valid in all orders of perturbation 
theory. This derivation establishes a relation between the parton densities 
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parametrized at n/ and rif + 1 light flavours. The consequences for the exist- 
ing parametrizations of the parton densities are discussed. Further we show 
that in charm electroproduction the exact and asymptotic expressions for 
the heavy-quark coefficient functions yield identical results for F 2)C (x, Q 2 , m 2 ) 
when Q 2 > 20 (GeV/c) 2 . We also study the differences between the FOPT 
and the VFNS descriptions for F 2iC (x, Q 2 , m 2 ). It turns out that the charm 
structure function in the VFNS is larger than the one obtained in FOPT 
over the whole Q 2 -range. Furthermore inspection of the perturbation series 
reveals that the higher order corrections in the VFNS are smaller than those 
present in FOPT for Q 2 > 10 (GeV/c) 2 . Therefore the VFNS gives a better 
prediction for the charm structure function at large Q 2 -values than FOPT. 



2 



1 Introduction 



The study of charm production in deep-inelastic electron-proton scattering 
and in photon-proton scattering provides us with important information 
about heavy-quark-production mechanisms. In particular, we can distin- 
guish between intrinsic- and extrinsic-charm production. In the latter case 
the charm quark only appears in the final state so the dominant subprocess 
is given by (virtual) photon-gluon fusion, which is the only reaction present 
on the Born level |rj . Hence one can measure the x-dependence of the gluon 
density G(x,[i 2 ), where x denotes the Bjorken scaling variable and /i stands 
for the factorization and renormalization scales. Next-to-leading (NLO) cor- 
rections 0, P|, to which also other processes contribute, reveal that this 
picture remains unaltered. 

Besides extrinsic-charm production one can also consider intrinsic-charm 
production 0], which is given by the flavour-excitation process. In this case 
the charm quark also appears in the initial state and it is considered to be a 
part of the hadronic wave function. Hence this quark is described by a parton 
density in the hadron, as in the cases of the other light flavours (u,d,s) and 
the gluon (g). Notice that in the context of perturbative QCD the flavour- 
excitation mechanism is always accompanied by the photon-gluon fusion pro- 
cess. However the latter can occur without the presence of intrinsic-charm 
production because one can always assume that the probability to find a 
charm quark inside the proton is equal to zero without violating the principles 
of perturbative QCD. The difference between both production mechanisms 
becomes very clear when one looks at differential distributions. In the case 
of extrinsic-charm production the photon-gluon-fusion subprocess predicts 
that both the charm quark and the charm anti-quark appear back-to-back in 
the final state. For intrinsic-charm production, where the flavour-excitation 
mechanism dominates, it turns out that either the charm quark or the charm 
anti-quark appears in the final state. Hence the transverse momentum of 
this heavy quark is not balanced by that of its antiparticle, in contrast to 
the situation in the photon-gluon fusion process. Recent experiments car- 
ried out at HERA at Q 2 = [| (photoproduction) or at 10 < Q 2 < 100 
(GeV/c) 2 || (electroproduction), where Q 2 denotes the virtual photon mass 
squared, give strong evidence for extrinsic-charm production. Therefore we 
only concentrate on the latter production mechanism in this paper. 

As has been mentioned above the extrinsic-charm mechanism can be stud- 
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ied in photoproduction (Q 2 = 0) and in electroproduction (Q 2 > 0). The 
photoproduction reaction has the advantage that the production rate is much 
larger than in the case of electroproduction because the latter rate is sup- 
pressed by the photon propagator, which decreases when Q 2 gets larger. 
In the context of perturbative QCD, however, the description of electro- 
production is easier due to the absence of the hadronic (resolved) photon 
component which contributes in photoproduction (for a discussion see @). 
Moreover electroproduction enables us to study the charm contribution to 
the deep-inelastic structure functions F2(x,Q 2 ) and Fl(x,Q 2 ). 

For the treatment of the charm component of the structure functions 
Fk, c (x,Q 2 ,ml) (k = 2,L) one has adopted two different prescriptions for 
extrinsic-charm production in the literature. The first one is advocated in 
|| where the charm quark is treated as a heavy quark and its contribution 
is given by fixed-order perturbation theory (FOPT). This involves the com- 
putation of the photon-gluon-fusion process mentioned above and its higher 
order corrections. The second approach is the so-called variable-flavour num- 
ber scheme (VFNS) ||. Here charm is treated in a similar way to a massless 
quark and its contribution is described by a parton density in the hadron, 
denoted by f^f™ s (x,fi 2 ). At first sight this looks similar to intrinsic-charm 
production. However there is one important difference. This becomes clear 
when one looks at eq.(10) of || where there exists a relation between jVfns 
and the gluon density originally appearing in the photon-gluon process. In 
this way there is an intimate relation between FOPT and the VFNS which im- 
poses some conditions on the charm-quark density. These conditions do not 
exist in the intrinsic-charm approach, where this density is just an arbitrary 
function fitted to experimental data. It is clear that the VFNS-procedure is 
not very well suited to describe charm production in the threshold region be- 
cause the photon-gluon fusion process requires that the charm component of 
the structure function vanishes above x = Q 2 / (Q 2 + 4 m 2 ), while this thresh- 
old condition gets lost in the x-behaviour of /,y FNS (a;, fJ* 2 )- On the other hand 
FOPT also has its drawbacks when Q 2 gets very large. The reason is that 
the heavy-quark coefficient functions calculated up to NLO in [[5J are domi- 
nated by large logarithms of the type \yl(Q 2 /m 2 ) hr 7 (/i 2 / m 2 .) when Q 2 3> m 2 



rpTOfl . Although it was shown in NLO [TT| that these logarithms lead to rather 
stable charm structure functions F i c (x,Q 2 ,m 2 ) with respect to variations in 
the factorization and renormalization scales /i, their size still warrants some 
special treatment (for a discussion of the scale dependence of the charm con- 
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tribution in the FOPT and VFNS approaches see also JT2[ , [0 ' ) • However 
as we will show in this paper the above criterion is not sufficient to decide 
about the rate of convergence of the perturbation series. It still may hap- 
pen that these large logarithms vitiate the perturbation series, in particular 
when corrections beyond NLO are included. Hence these logarithms have to 
be resummed using the standard techniques of the renormalization group. 

This resummation proceeds as follows. First one has to add the light- 
parton component of the structure functions Fk(x, Q 2 ), k = 2, L to the charm 
contributions Fk }C {%, Q 2 , m 2 ) both presented with three light flavours. The 
light component consists of the light-parton (u,d,s,g) densities convoluted 
with the light-parton coefficient functions. The same densities are used for 
the charm component where they are convoluted with the heavy-quark coeffi- 
cient functions. Then one takes the limit Q 2 — > oo and performs mass factor- 
ization in order to remove the above mass singular logarithms in m c from the 
asymptotic heavy-quark coefficient functions. The latter can be written as a 
convolution of the heavy-quark operator-matrix elements (OME's), contain- 
ing these mass-singular logarithms, and the light-parton coefficient functions 
which are finite in the limit m c — > 0. The light-parton densities are then 
modified by multiplying the original u,d,s and g densities by the heavy-quark 
OME's, where one has to introduce a new parton density to represent the 
charm quark. Notice that the latter is not an arbitrary function since, via the 
heavy-quark OME's, it depends on the original densities represented by the 
three light flavours (u,d,s) and the gluon (g) in the three-flavour scheme. The 
whole procedure leads to the structure functions Fi and F ifi which are now 
represented in the four-flavour scheme. The latter are expressed into four 
light-flavour densities (including charm) and the gluon density convoluted 
with the light-parton coefficient functions. Further the heavy-quark coeffi- 
cient functions due to the charm have completely disappeared. Although 
these formulae are similar to those obtained for intrinsic-charm production 
their origin is completely different. Here we want to emphasize again that 
the VFNS is derived from FOPT, where, in the large-Q 2 limit, one has re- 
summed the large logarithmic terms in the heavy-quark coefficient functions 
in all orders of perturbation theory. This procedure, which is only described 
up to leading order (LO) in ||, will be generalized to all orders in QCD 
perturbation theory in Section 2. 

In Section 3 we will study at which Q 2 the large logarithms in the heavy- 
quark coefficient functions actually dominate the charm component of the 
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structure functions. We will investigate the x- and /^-dependence of the 
charm density in the VFNS approach and compare it with the existing charm 
densities in the literature. We will also study the differences between the 
FOPT approach and the VFNS approach for the charm component of the 
deep-inelastic structure functions in the x- and Q 2 -range explored by present 
experiments. Further we investigate which approach is more stable with 
respect to higher order corrections to the charm structure functions. In 
Appendix A we present some heavy-quark OME's which were not previously 
calculated in the literature. In Appendix B we list all renormalized heavy- 
quark OME's which are needed for our analysis in Section 3. 
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2 Derivation of the VFNS representation of 
the structure functions 



In this section we present the variable-flavour number scheme (VFNS) rep- 
resentation of the structure functions F^(x, Q 2 ) in all orders of perturbation 
theory. Our results hold for any species of heavy quark although at present 
collider energies the VFNS is only interesting for the charm quark. As men- 
tioned in the previous section logarithms of the type \n l (Q 2 /m 2 ) In 7 (/i 2 /m 2 ) 
arise in the heavy-quark coefficient functions when Q 2 ^> m 2 and we work 
in fixed-order perturbation theory (FOPT). Here m stands for the mass of 
the heavy quark, denoted by H in the subsequent part of this section. When 
going from the FOPT representation for the structure functions to that of 
the VFNS one has to remove these mass-singular logarithms from the heavy- 
quark coefficient functions. This is achieved using the technique of mass 
factorization, which is a generalization to all orders in perturbation theory 
of the procedure carried out up to lowest order (LO) in Sections II D and II 
E in || . Although this procedure resembles the usual mass factorization of 
the collinear singularities, which appear in the partonic structure functions 
or partonic cross sections, it is actually much more complicated. This is due 
to the presence of the light (u,d,s and g) partons, as well as of the heavy (c,b 
and t) quarks, in the Feynman diagrams describing heavy-flavour production. 

In the calculations of the heavy-flavour cross sections the light partons 
are usually taken to be massless, whereas the heavy quarks get a mass m^O. 
Note that the mass is defined by on-mass-shell renormalization. When Q 2 ^> 
m 2 two types of collinear singularities appear in the partonic contributions. 
One type can be attributed to the light partons. In this case the collinear 
divergences can be regularized using various techniques. The most well- 
known among them is n-dimensional regularization. The second type can 
be traced back to the heavy quark and the singularity manifests itself as 
m — > in the large logarithmic terms mentioned above. Beyond order a s 
both types of singularities appear in the partonic cross sections and in the 
heavy-quark operator matrix elements (OME's). The latter show up in the 
mass-factorization formulae, and they are needed to remove the mass-singular 
terms (as m — > 0) from the heavy-quark coefficient functions. Hence the 
mass factorization becomes much more complicated than when we only have 
to deal with collinear divergences due to massless partons. Therefore we first 
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present the mass factorization with respect to the latter before we apply this 
technique to the mass singularities related to the heavy-quark mass m. 

Consider first deep-inelastic electron-proton scattering in which only light- 
partons show up in the calculation of the QCD corrections. The deep- 
inelastic structure function Fj(n/,Q 2 ) can be expressed as follows 



F t (n f ,Q 2 ) = - E 
n f fc=i 



/ 

2 



±{n f ) ® jfjn f , ^, /i 2 ) + G(n f ) ® ^ g (n f , ^, /i 2 ) 



+n / A fc (n / )®^ s (n / ,^, / u 2 ) 



(2.1) 



In this equation the charge of the light quark is represented by e& and rif 
denotes the number of light flavours. Further <g> denotes the convolution 
symbol and, for convenience, the dependence of all the above quantities on 
the hadronic scaling variable x and the partonic scaling variable z is sup- 
pressed. The bare quantities in (2.1) are indicated by a hat in order to 
distinguish them from their finite analogues, which emerge after mass fac- 
torization. Starting with the parton densities, X( n /) an d G{nf) denote the 
singlet combination of light-quark densities and the gluon density respec- 
tively, with rif light flavours. The former is given by 



E(n/) = E [fiM + fiM] . (2-2) 
i=i 



where fiijtf) and fj{rif) stand for the light-quark and light-anti-quark den- 
sities respectively. The non-singlet combination of light-quark densities is 
given by 

A fc (n/) = f k (n f ) + f- k (n f ) - ±-±(n f ) . (2.3) 

The QCD radiative corrections due to the (virtual) photon light-parton 
subprocesses are described by the partonic structure functions J-ij(nf) (i = 
2, L; I — q, g) where I stands for the parton which appears in the initial state. 
As with the parton densities they can also be classified into singlet, non- 
singlet and gluonic parts. Furthermore we assume that coupling-constant 
renormalization has been performed on the partonic (bare) structure func- 
tions, which is indicated by their dependence on the renormalization scale 
fi. However they still contain collinear divergences which, for convenience, 
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are regularized by taking the external momentum p of the incoming parton 
off- mass-shell (p 2 < 0). Notice that these divergences do not show up in the 
final state because the deep-inelastic structure functions are totally-inclusive 
quantities. 

The reason we choose off-shell regularization is that it allows us to dis- 
tinguish between the collinear divergences coming from the massless partons 
and from the heavy quarks. The former divergences are the ln l (— p 2 /p 2 ) 
terms in the perturbative expansion, whereas the latter are the ln*(/z 2 /m 2 ) 
terms. 

We can now also reexpress Eq. (2.1) in finite quantities so that the 
collinear divergences are absent. This is achieved via mass factorization 
which proceeds as follows 
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and 



*ff(»/. = < S (>V. g) • <$?(»/■ %) • (2.4) 

In the above expressions (^(i = 2, L; k = q,g) denote the light-parton co- 
efficient functions and the represent the renormalized operator-matrix 
elements (OME's) which are defined by 

M n /> ^) =< Hp)\Oi(0)\Hp) >, (l,k = q,g). (2.6) 

Here 0/ are the renormalized operators which appear in the operator-product 
expansion of two electromagnetic currents near the light cone. The product 
of these two currents is sandwiched between proton states and its Fourier 
transform into momentum space defines the structure functions in (2.1). Like 
the other quantities and A ik can be divided into singlet and non-singlet 
parts. The scale p appearing in Eqs. (2.4)-(2.6) originates from operator 
renormalization as well as from coupling-constant renormalization. Notice 
that the operator-renormalization scale is identical to the mass-factorization 
scale. Using the mass-factorization relations in Eqs. (2. 4), (2. 5) we can cast 
the hadronic structure functions Fi(nf,Q 2 ) (2.1) into the form 



i n f 

Fi(n f ,Q 2 ) = -E4 
n f k=i 



nn f ^ 2 ) ® Cf tq (n f , ^) + G(n f ,p 2 ) ® ^) 
p P 
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+n f A k (n f ,^)^C^(n f ,%)}, (2.7) 

where the finite (renormalized) parton densities S, A and G are expressed 
in the bare ones S, A and G in the following way 

A k (n f ,^ 2 ) = Al s (n f ,^) ® A k (n f ) , (2.8) 
/i 2 ) = ^(n/, ^) ® S(n/) + 4,(n/, ^) ® G(n,) , (2.9) 



and 



G(n/, ^ 2 ) = ^(n,, ® t(n f ) + (n,, ® G(n f ) . (2.10) 

All quantities given above satisfy renormalization group equations (RGE's). 
Here we are only interested in those RGE's for the OME's and the parton 
densities. Define the differential operator D as 

D = fi— + P(n f ,g)—, g = g(n f ,fi 2 ), (2.11) 
where a s = g 2 /(4ii), then the OME's satisfy the following RGE's 

DAl S (n f , = - 7 -(n /} ® A^(n f , , (2.12) 

S = - E 7&(n/) ® A%(n f , . (2.13) 

" k=q,g " 

Here 7 y denote the anomalous dimensions corresponding to the operators 
Oi. They can be expanded as a perturbation series in a s . In Bjorken x- 
space there exists a relation between the anomalous dimensions 7^ and the 
DGLAP splitting functions, denoted by P^, which is given by 

7i >/) = -P«(n/) . (2.14) 

Notice that this relation only holds for twist-two operators. From the above 
equation one infers that the 7^, which are the residues of the ultraviolet 
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divergences in the unrenormalized OME's, have just the opposite signs to 
those of the P^. The latter show up in the partonic quantities ^Fi,i(nf) in 
(2.1) and they represent the residues of the collinear divergences. We will 
return to this relation (2.14) when we discuss the heavy-quark OME's. The 
finite (renormalized) parton densities satisfy the RGE's 

DA k (n f ,fi 2 ) =-7£ S (n/)® A k (n f ,fj, 2 ), (2.15) 

DH{n f ,n 2 ) = -ll{n f ) ® S(n /)/U 2 ) - ^{nj) ® G(n f ,fi 2 ) , (2.16) 

and 

DG(n f , n 2 ) = -i s gq (n f ) ® S(n /; fi 2 ) - >f gg {n f ) ® G(n f , fi 2 ) . (2.17) 

From these equations we can derive the Altarelli-Parisi equations. 

Before we add the heavy-quark contributions to the deep-inelastic struc- 
ture functions (2.7) it is convenient to split the singlet quantities jFf q1 Cf q 
and A qq into non-singlet and purely-singlet parts, namely 

*&=3y+3?, (2-i8) 

<t = C + C£, (2-19) 

and 

A qq = Aqq + A qq . (2.20) 

This decomposition facilitates the mass factorization of the heavy-quark coef- 
ficient functions and can be explained as follows. If we calculate the diagrams 
contributing to the parton subprocesses with a quark in the initial state, the 
resulting expressions have to be projected on the singlet and non-singlet chan- 
nels with respect to the flavour group. The latter projection leads to ^ N g s . 
However the singlet part i.e. Tf can be split into two types of contribu- 
tions. The first one is equal to Pf q whereas the second one is represented by 
fif q - The purely-singlet partonic structure function arises from the Feynman 
graphs where the projection on the non-singlet channel yields zero, so that 
only singlet contributions remain. They are characterized by those Feynman 
graphs in which only gluons are exchanged in the t-channel. Such graphs 
show up for the first time in two-loop order. An example is given in Fig. 
1. The same characteristics also hold for (see Fig. 2) and the resulting 
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coefficient function Cf q . Another important feature is that the purely-singlet 
quantities are proportional to the number of light flavours rif. This property 
is shared by the gluonic quantities fif , Cf and A q . The proportionality to 
rif can be traced back to the fact that in the case of Ffq, ^Ff :9 the virtual 
photon is attached to the light-quark loop (see e.g. Fig. 1) and one has to 
sum over all light flavours. In the case of A^ q and A qg this is due to the 
insertion of the operator vertex into the light-quark loop where also a sum 
over all light flavours has to be carried out. Because of the mass-factorization 
relation (2.5) the proportionality to rif is transferred to the coefficient func- 
tions Cf t g, Cf g and the anomalous dimensions 7^, 7^ . To facilitate the mass 
factorization of the heavy-quark coefficient functions it is very convenient to 
extract this overall factor of rif from the quantities above so we define 

T ijk = n f T iyk , (2.21) 

where 

rp *rPS pYS . rp ^-S /rS . /o 99^ 

1 i,q r i,g j ^i, q j 1 i,g ~ S \,gi ^i,g j l Z ' ZZ J 

and 

Rij = n fRij , (2.23) 

where 

R qq = A™ n l s ; R qg = A s qg n s qg . (2.24) 

Besides the overall dependence on rif, which we have extracted from the 
quantities defined by T^k and R^ above, there still remains a residual depen- 
dence on rif in T i>k and R^. The latter dependence originates from internal 
light-flavour loops which are neither attached to the virtual photon nor to the 
operator- vertex insertions. Since there are more contributions to T itk and R^ 
which are not due to these light flavour-loops it is impossible to extract an 
overall factor rif anymore from the quantities indicated by a tilde. Therefore 
the dependence of the latter on rif means that it can be only attributed to 
internal light-flavour loops. Using the above definitions one can now rewrite 
Eqs. (2.1) and (2.7) and the results become 



n f 

F t (n f ,Q 2 ) = Y.e\ 
k=i 



~PS , Q2 ^ ~ s Q2 

Z(n f ) <g> T iA [rif, — , Li 2 ) + G(n f ) ® ^{rif, -j, fi 2 ) 



Q 2 



+{/*("/) + h{nf)}®^{rif^^ 2 ) 



(2.25) 



12 



F t {n f) Q 2 ) = E4 nn f ,^)0^(n f ,^)+G(n f ,^)^ef :g (n f ^) 
k=i L Z 1 Z 1 

+ {/ fe (n /!/U 2 ) + A(n / ,/i 2 )}®C^ s (n / ,^)j , (2.26) 



with the relation 



/,(%,/i 2 ) + A(%,/i 2 ) = <(n/,^)® {/*("/) + A("/)} 

+ A ^{^ S ® Z(n f ) + ^(n,, ® G(n,) . 



Using Eqs. (2. 5), (2. 9) and (2.10) one can now write the mass factorization 
relations for the quantities indicated by T itk and Rij in (2.21)-(2.24). The 
same can be done for the RGE's which can be derived from (2. 13), (2. 16) and 
(2.17). Since this derivation is easy it is left to the reader. Here we only 
want to report the RGE for the left- hand- side of (2.27) which follows from 
(2.15) and (2.16). It is given by 



D[f k (n f ^ 2 )+f k (n f ^ 2 )] = ^ 8 (n,)®[f k (n f ,iJL 2 ) + f- k (n f ,iM 2 ) 



After having presented the formulae needed for the mass factorization of 
the light-parton structure functions jF ifl with the corresponding RGE's we 
want to deal in a similar way with the asymptotic heavy-quark coefficient 
functions where the large logarithmic terms depending on the heavy-quark 
mass m have to be removed. For that purpose we have to add the heavy- 
quark contribution to the deep-inelastic structure function Fi(rif,Q 2 ) (2.26) 
which is equal to 



(2.27) 



% P <>/) ® /i 2 ) - f qg (n f ) <g> G{n f , /i 2 ) . 

(2.28) 



F itH (n f ,Q 2 ,m 2 ) = 



n f r C) 2 m 2 

•£4 E{n„ ,*>) « ™) 

k=i L P 




+ {h(n f ^ 2 ) + h(n f ,^ 2 )} ® L™(n f , 
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+e 2 H 



(2.29) 



where e# stands for the charge of the heavy quark denoted by H. Further 
L i>k and H i k {i — 2, L; k — q,g) represent the heavy-quark coefficient func- 
tions. Like in the case of the light-parton coefficient functions d$ they can be 
split into singlet and non-singlet parts. The former can be decomposed into 
non-singlet and purely-singlet pieces in a similar way as given in (2.19). The 
distinction between and can be traced back to the different (virtual) 
photon-parton heavy-quark production mechanisms from which they origi- 
nate. The functions are attributed to the reactions where the virtual 
photon couples to the light quarks (u, d, and s), whereas the describe 
the interactions between the virtual photon and the heavy quark. Hence 
and k in (2.29) are multiplied by ef and e 2 H respectively. Moreover, when 
the reaction where the photon couples to the heavy quark contains a light 
quark in the initial state, then it can only proceed via the exchange of a 
gluon in the t-channel (see Fig. 1). Therefore H ijq is purely-singlet and a 
non-singlet contribution does not exist. This is in contrast with Lj i? which 
has both purely-singlet and non-singlet contributions. 

For Q 2 ^> m 2 the heavy-quark coefficient functions take their asymp- 
totic forms, which are dominated by the large logarithms mentioned at the 
beginning of this section. If these terms become too large they vitiate the 
perturbation series so that a resummation via the RGE is necessary. Before 
this resummation can be carried out we have first to perform mass factor- 
ization to remove the mass-singular terms ln*(/x 2 /m 2 ) from the asymptotic 
heavy-quark coefficient functions. This can be done in a similar way as 
shown for the partonic structure functions in Eqs. (2. 4), (2. 5), where now 
ln l (— fi 2 /p 2 ) are replaced by ln l (/i 2 /m 2 ). In the VFNS approach we require 
that the following relation holds 

Fi(n f , Q 2 ) + lim \Fi, H (n f , Q 2 , m 2 )} = F^ FNS (n f + 1, Q 2 ) . (2.30) 

The above formula implies that, at large Q 2 , the left-hand-side can be writ- 
ten as a structure function containing densities and coefficient functions cor- 
responding to only light partons, in which the number of light flavours is 
enhanced by one. In Fj VFNS (n/ + l,Q 2 ) the light-parton densities are mod- 
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ified with respect to those appearing in Fi(n/,Q 2 ) in (2.26) and the heavy- 
flavour component F it H(Q 2 ,m 2 ) in (2.29) because the former have absorbed 
the logarithms In* (/i 2 /m 2 ) coming from the heavy-quark coefficient functions. 
Moreover a new parton density appears in i^ VFNS corresponding to the heavy 
quark H which is now treated light quark. 

If relation (2.30) holds one has to impose the following mass-factorization 
relations for the heavy-quark coefficient functions. In the case of the coeffi- 
cient functions L ijk we get for Q 2 3> m 2 

C™(n f ) + L™(n f ) = Al s H (n f ) ® C™(n f + 1) , (2.31) 



<?>/) + ^>/) = + n S A 2n(n f ) + A%(n f )\ ® C™(n f + 1) 



iNS 



PS 



IPS, 



TPS/ 



+ A2„(n f ) ® CjJ(n/ + 1) + ^ iH (n/) ® C?>/ + 1) , 

(2.32) 



and 



<?>/) • f ■).,("/) = ^(n/) ® C(n/ + 1) + ^(n/) ® Cf>/ + 1) 



+ 



+ I ® + 1) • (2.33) 



For the heavy quark coefficient functions we obtain for Q 2 ^> m 2 



H i,q( n f) = A™(n f )®\C™(n f + l)+C™(n f + l) 



I PS 



NS/ 



5T-S/ 



+ 



+^Lh(«/)®C?>/ + 1)> 



(2.34) 



and 



+A» Bg (n f ) ® lC(n f + 1) + Cf>/ + 1) 



PS/ 



(2.35) 



In the above equations A Hk denotes the heavy-quark OME defined by 



A Hk (n f ,^) =< k(p)\O H (0)\k(p)>, 



(2.36) 
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which is the analogue of the light-quark OME's in (2.6). The quantities 
Aik,H(nf,fi 2 /m 2 ), which also appear above, represent the heavy-quark-loop 
contributions to the light-quark and gluon OME's defined in (2.6). 

Although there exists some similarity between the mass factorization of 
m (2-4), (2. 5) and the ones given for the heavy-quark coefficient functions 
above we also observe a striking difference which makes the proof of Eqs. 
(2.31)-(2.35) and consequently of Eq. (2.30) much harder. The main differ- 
ence is that the light partons appear in the initial as well as in the final state 
of the subprocesses contributing to Ti t i whereas the heavy quark only shows 
up in the final state of the reactions leading to L^H^i (2.29). This implies 
that there is no analogue for !F^h in Eq. (2.29) meaning that the coefficient 
functions L^h and H^h do not appear in the latter equation. The same holds 
for the bare heavy-quark density /# which has no counter part in (2.26) and 
(2.27) either. Hence the mass factorization relations in (2.31)-(2.35) are much 
more cumbersome than those presented in (2. 4), (2. 5). Therefore the proofs 
of the former relations and implicitly of Eq. (2.30) become more involved. 
We have explicitly checked that the above relations hold up to order a 2 using 
the asymptotic heavy-quark coefficient functions in [10] and the heavy quark 
OME's in Appendix B. 

Yet another complication arises when we consider the OME's A^k and 
Aik,H defined above. Contrary to the light-parton OME's in (2.6) which 
only depend, apart from fi 2 , on one mass scale p 2 the former also depend 
on the mass scale m 2 due to the presence of the heavy quark. Therefore 
the unrenormalized expressions of Ahu and A^^h contain besides ultraviolet 
(UV) divergences two types of collinear divergences which are represented 
by ln*(/x 2 /p 2 ) (light partons) and ln*(/x 2 /m 2 ) (heavy quark) respectively. The 
singularity at p 2 = shows up because external massless lines represented by 
k = q, g in (2.36) are coupled to internal massless quanta. This phenomenon 
shows up for the first time in order a 2 see ||10[ . Notice that the singularities 
at p 2 = also appear in the partonic quantities leading to the heavy-quark 
coefficient functions before they are removed by mass factorization as outlined 
in the beginning of this section. Therefore we also have to subtract the 
collinear singularities in p 2 = from the unrenormalized OME's Ank,Aik,H 
in addition to the UV divergences. This twofold subtraction leads to two 
different scales in the renormalized OME's in (2.31)-(2.36) which are called 
the operator-renormalization scale and mass-factorization scale respectively. 



16 



Usually these two scales are set to be equal and they are denoted by one 
parameter /i 2 . However the appearance of these two scales results in more 
complicated RGE's for the heavy-quark OME's in comparison with those 
presented for the light-parton OME's in (2. 12), (2. 13) as we will see below. 
In the actual calculations of A Hk ,Ai k ^ H in [|l(J we have put p 2 = in (2.36) 
because we did the same for the partonic quantities computed in leading 
to the heavy-quark coefficient functions L^ j, H it i. Hence we had to adopt 
n-dimensional regularization for the UV as well as the collinear singularities 
in p 2 = 0. In this way both are represented by pole terms of the type 
l/(n — 4) J . If the latter are removed in the MS-scheme the OME's Auk and 
Aik,H automatically depend on one scale \i. 

Another comment we want to make is that a s appearing in the heavy- 
quark coefficient functions is renormalized at rif flavours. This means that 
the renormalization of the coupling constant is carried out in such a way that 
all quarks equal to H and heavier than H decouple in the quark loops con- 
tributing to aig. In the right- hand- side of Eqs. (2.31)-(2.35) the decoupling 
of the heavy quark H has been undone so that here a s depends on rif + 1 
flavours. This decoupling gives rise to additional logarithms of the type 
ln l (/z 2 /m 2 ), which are cancelled by the OME's of the type A^h- However 
explicit indication of this procedure further complicates our mass factoriza- 
tion relations above, which we want to avoid, so we have to bear in mind 
that this decoupling is implicitly understood. 

If we now substitute Eqs. (2.31)-(2.35) into the left-hand-side of (2.29) 
and rearrange some terms then we obtain the expression for F fc VFNS as pre- 
sented in (2.26), where rif — > n/ + 1. This implies that the new parton 
densities taken at rif + 1 light flavours can be expressed in terms of those 
given for rif flavours. The original rif light-flavour densities get modified so 
that for k — 1, . . . , rif 

n 2 

|NS / h 1 



f k (n f + l,fi 2 ) + f- k (n f + l,n 2 ) = A™ H (n f , ^) ® [f k (n f ,fi 2 ) + f- k (n f ,^ 

+ A ™H(rif,^)®X(rif, f i 2 ) 
+ A l,H(n f ^)®G{n f ^ 2 ) ) (2.37) 
whereas the parton density of the heavy quark can be expressed in the original 
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light flavours in the following way 

fH+it(nf + 1, V 2 ) = fn f +i(n f + 1, /i 2 ) + f^j+i(n f + 1, fi 2 ) 

2 2 

^K( n /> ^^Z"" 2 ) + A>( n /> £2) ® G ( n f^ 2 ) ■ 

lit lit 

(2.38) 



Comparing the above expression with Eq. (2.27) we observe that the first 
term in (2.27) has no counter part in (2.38). This is because we have no 
bare heavy-quark density unless one assumes that there already exists an 
intrinsic heavy-quark component of the proton wave function. The singlet 
combination of the quark densities becomes 



S(n / + 1^ 2 )= ]T [f k ( nf + l,fi 2 ) + f- k (n f + l,fi 2 ) 
k=i 

2 22 
= A TA n f> £2) + M£h("/> £2) + ^{rif, ^) 
2 

" A* ( n /> £2) + ^ s ( n /' £2") 



+ 



S(n/,/i 2 ) 
(2.39) 



The non-singlet combination A fc (n/ + 1) is defined in an analogous way as 
in (2.3) and it reads for k — 1, . . . , rif + 1 



A k (n f + 1, /i 2 ) = / fc (n, + 1, /i 2 ) + /^n, + 1, /1 2 ) — + 1, fi 

n f + 1 



(2.40) 



Finally the gluon density for rif + 1 light flavours is 

G(n/ + l,/x 2 ) = ^(n;,^ 2 )®^^/,^ 2 ) 

+Af-„(n / ,Ai 2 )®G(n/,Ai 2 )- 



(2.41) 



The old as well as the new parton densities have to satisfy the momentum 
sum rule 



Jo 



dxx 



E(rif,x, /i 2 ) + G(n f ,x, /i 2 ) =1, 



(2.42) 
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for any rif. This implies that the OME's Auk, Au,h have to satisfy two 
relations 



2 2 

I dxx A™ H (n f ,x,^)+n f A™ H (n f ,x,^) 



2 2 



and 



I dxx n f A l 9 ,u(nf ) x ) ^) + A s Hg (n f ,x, 



(2.43) 



(2.44) 



which can be checked up to second order using the results in Appendix B. 
The OME's A Hk and A klH satisfy the following RGE's 



DA PS 



(7,T + -/7 PS 



IPS i „,s 



4 s 



Iqq 1 '"J iqq 

-(1h% + Thh) ® A ™ - f Hg ® A b gq:H 
-f H S®{Al S H + n f A™ H ), 



(2.45) 



UA qq,H 



lqq,H 



aNS 
A qq,H » 



4 s 
Hg 



(2.46) 
(2.47) 



D4 PS 

U/1 qq,H 



7 S 



A qg,H \lqq,H + n flqq,H 



PS 



[PS 



5 

4 s 



(2.48) 



D4 S 



= fqg ® (^H + n/ijy + ^ ^ 
- (tJ^J + "/T^S + + n flqq) ® 



7 S 

'9S 



...... 



;PS ^ 4S 



(2.49) 
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DA s gqiH = (^ q s + na™)®A s gqiH + 1 s gq ®A s ggiH 

" bl,H + ll) ® - 7 9 S h ® A™ , (2.50) 



DA Lh = n f f qg ® A s gq , H - ^ H ® A% iH 

S i _,S i 4S „,S 4S 



and 



-«/ {r 9q ,H + r gq ) ® - %h ® A b Hg . (2.51) 

The above RGE's are much more complicated than those written for 
A k i in Eqs. (2. 12), (2. 13). This is due to the fact already mentioned be- 
low (2.36) that fj, represents the operator-renormalization scale as well as 
the mass-factorization scale. In Eqs. (2. 12), (2. 13) it only stands for the 
operator-renormalization scale. The anomalous dimensions coming from op- 
erator renormalization carry a minus sign whereas those associated with mass 
factorization have a plus sign in front of them. This difference in sign can be 
traced back to Eq. (2.14) where it was stated that the residues of the ultravio- 
let divergences are just the opposite of the ones corresponding to the collinear 
divergences. Furthermore one has to bear in mind that the residues of the 
ultraviolet and collinear divergences are equal to the anomalous dimensions 
coming from operator renormalization and mass factorization respectively. In 
Eqs. (2.45)-(2.51) the anomalous dimensions 7^ have a plus sign on account 
of the collinear divergences occurring in A k i, where the partons indicated by 
k, I are massless. However the anomalous dimensions 'Jhh, Thz, Iih and ^ki,H 
carry a minus sign because the mass of the heavy quark m prevents these 
OME's from being collinearly divergent so that we only have to deal with 
ultraviolet singularities. 

Using the above equations and (2.15)-(2.17) one can derive the RGE's for 
the new parton densities appearing in _pY FNS (2.30). For k — 1, . . . , rif + 1, 
including the heavy-quark flavour, the RGE reads 

D 



f k (n f + l,n 2 ) + f- k (n f + l,Li 2 ) 
-7j?("/ + 1) ® [/*("/ + 1, V 2 ) + h(n f + 1, /i 2 ) 



-^ q s (n f + 1) (g) Y.{n f + l,/i 2 ) - f qg {n f + 1) <g> G(n f + l,/i 2 ) .(2.52) 
The non-singlet combination of the quark densities satisfies the RGE 

DA k (n f + 1, /1 2 ) = - 7g N g S (n/ + 1) ® A k (n f + 1, /1 2 ) . (2.53) 
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The singlet combination of the quark densities satisfies the RGE 



DE{n f + l,fi 2 ) = -7^(71/ + 1) <g> £(7i/ + 1,// 2 ) 

-7 9 S 9 (^ + l)®G(n / + l,/i 2 ), (2.54) 

and the gluon density is given by 

DG(n f + l,ti 2 ) = -7 g s ? (n/ + l)®£(n/ + l,/x 2 ) 

-7 s S 9 (n/ + 1) ® G(n f + 1, /. 2 ) . (2.55) 

In the above equations we have used the identities 

7y(»v) + lij,H(n f ) = >y tj (n f + 1) , (2.56) 

lW H (n f + 1) = 7, N >/ + 1) J 7$r(»/ + 1) = % P >/ + 1) , (2-57) 
% P K% + 1) = 7^(n, + 1) = 7 P >/ + 1) , (2.58) 

l S gH (n f + 1) = 7 9 S >/ + 1) ! f Hg (n f + 1) = % s 9 (n/ + 1) , (2.59) 

because the anomalous dimensions do not depend on the mass m of the 
heavy-quark H. A comparison of the above RGE's with those presented for 
the light-parton densities in Eqs. (2.15)-(2.17) reveals that they are exactly 
the same in spite of the fact that there is no counterpart of the bare heavy- 
quark density in the derivation of Eqs. (2.52)-(2.55). 

All perturbative quantities which appear in this section are now available 
up to 0(ot 2 s ). This holds for the anomalous dimensions 7m JT5], the massless- 



parton coefficient functions Ci t k |J6 and the heavy-quark coefficient functions 
L it k, H itk [§. The OME's A^ q , A s Hg and A^ H are computed up to 0(a 2 s ) 
and listed in unrenormalized form in Appendix C of [[TIJ. We still require 
H (Fig. 3) and Ag H (Fig. 4) which are calculated up to the same order 
in this paper. Exact expressions for the unrenormalized OME's can be found 
in Appendix A. Notice that up to second order in a s both A qg H and A^ H 
are zero. The renormalized (finite) expressions for all these OME's, which 
we will use in the next section, are presented in Appendix B. 
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After having found the representation of the heavy-quark density fn+R 
in (2.38) and the RGE in (2.52) which determines its scale evolution we can 
write down the charm component of the deep-inelastic structure function in 
the VFNS representation. The latter is given by 



F^ S (n f + l,Q 2 



W(n/ + l,/i 2 )®C(n / + l,^) 



+S(n / + 1,^)0^(^ + 1,^ 



(2.60) 



The expression above satisfies DF^™ S = (see (2.11)), so it is renormaliza- 
tion group invariant, which means that it is scheme independent and becomes 
a physical quantity. Notice that even though the form of F^™ s is the same as 
the one presented for intrinsic heavy-quark production, their origins are com- 
pletely different. The VFNS result, Eq. (2.60) is derived from FOPT when 
the exact heavy-quark coefficient functions are replaced by their asymptotic 
expressions taken at Q 2 3> m 2 (see (2.30)). Then we have performed mass 
factorization and absorbed the mass singularities with respect to the heavy- 
quark mass m into the light-parton densities taken at the original rif light 
flavours (see Eqs. (2.37)-(2.41)). This procedure leads to relations between 
the parton densities in the n/ and rif + 1 flavour schemes. In particular this 
applies to the heavy-quark density fn+B i n (2.38). For intrinsic heavy-quark 
production the origin of Eq. (2.60) is completely different. In this case it is 
not derived from perturbation theory and therefore there does not exist any 
relation between the heavy-quark density and the light-parton densities like 
in Eq. (2.38). Actually the intrinsic-charm density is determined by a simple 
fit to experimental data. 

Furthermore we want to emphasize that the VFNS approach is only valid 
for totally inclusive quantities like structure functions since the logarith- 
mic terms ln l (Q 2 /m 2 ) ln J (n 2 /m 2 ) only appear in the asymptotic form of the 
heavy-quark coefficient functions in this case. Hence expression (2.60) is just 
an alternative description for FOPT when the production of the heavy quarks 
occurs far above threshold where Q 2 S> m 2 . The only difference between the 
FOPT and the VFNS descriptions is that the large logarithms have been re- 
summed in the latter approach so that one gets an improved expression with 
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respect to normal perturbation theory in the large Q 2 -region. However on the 
level of differential distributions the large logarithms of the kind given above 
do not show up in the perturbation series. Therefore in this case the VFNS 
approach cannot be applied. As has been mentioned in the introduction one 
can only distinguish between intrinsic and extrinsic heavy-quark production. 

Finally we want to comment about the work in || where one has proposed 
the idea of the VFNS approach. In particular we want to make some remarks 
about equation (9) in ||(ACOT) which is similar to our equation (2.60). 
Using the notations in the latter reference this equation reads 

E = fS ® E "bq - ® f? {l) ® E <42 } + f 9 N ®T, <4? } , 

A A A A 

(2.61) 

where we have summed over all helicities of the virtual photon denoted by 
A. Further we have corrected a misprint because /^ (0) in (9) should read 
To translate (2.61) into our language we have to make the following 
replacements 

B^7*; Q^H- E^iv-i^° T , (2-62) 

A 

E 42 } - 4 C (0) = -*); E < - e%H?jv , (2 . 63 ) 

A A 



fQO) 



A 



s,(i) 

Hg J 



(2.64) 



and 



J N 



G. 



/iV — *■ fu+H j 

Hence Eq. (2.61) can be written in our notation as 



pACOT _ 2 



w®C (0) +G®a s,(1) 



+4 



G ® {fzg. (1) - C. 



5S,(1) 



.4 



S,(l) 



NS,(0) 



V.(/ 



(2.65) 



(2.66) 



Comparing this expression with ours in (2.60) we notice the following differ- 
ences. First Eq. (2.61) and therefore Eq. (2.66) in |J was only derived in 
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lowest order (LO), whereas Eq. (2.60) is valid in all orders of perturbation 
theory. Second Eqs. (2.61) and (2.66) are not valid in next-to-leading order 
(NLO) but this can be repaired by replacing Cf^ by Cf^ +(a a /47r) Cfj? 
in (2.66). However the most important difference is that the last term be- 
tween the square brackets in (2.66) has no counter part in our expression 
(2.60). This is because we took the limit Q 2 — > oo and dropped all terms 
proportional to (m 2 /Q 2 ) 1 in the asymptotic heavy-quark coefficient functions. 
This is revealed by the mass-factorization relation (2.35) for Hf , which only 

holds in the limit Q 2 ^> m 2 . However the exact expressions for , 

and A*gg can be found in the literature (see ||, fl6| and £L0|). Expanding 
the latter in powers of m 2 /Q 2 we obtain 



„S,(1) /5S,(1) IS,(1) 
n i,g u i,g ^Hg 



Oi 



00 m 2 . 



1Q 2 



c 



NS,(0) 



l,q 



.(0 



ln( 



Q 



m 



(0 



(2.67) 



When Q 2 — > 00 (2.67) vanishes, as was already expected from (2.35). The 
motivation for the above expression, included in [[|, was to get a better 
stability of F^™ s in the threshold region with respect to variations in the 
factorization scale. However we have dropped these type of contributions as 
shown in Eq. (2.67) in our representation for F^™ s (2.60) for theoretical as 
well as practical reasons. The theoretical argument is that mass factorization 
does not apply to terms proportional to (m 2 /Q 2 ) 1 for / > 1 so that Eq. (2.66) 
or (9) in |J cannot be generalized to higher orders in a s . The practical reason 
is that the inclusion of higher order corrections automatically improves the 
stability of the perturbation series with respect to variations with respect to 
the factorization scale as is shown in [11], [12|], |13[. 



Summarizing this section we have presented expressions for the total 
structure functions Fj VFNS (?2/ + 1,Q 2 ) (2.30) as well as for its heavy-quark 
component F^ NS (nj + 1,Q 2 ) (2.60) in the context of the variable-flavour- 
number scheme (VFNS), which are valid in all orders of perturbation theory. 
Starting from the structure functions Fi(rif,Q 2 ), expressed in the light-parton 
densities and coefficient functions, we could show that by adding the heavy- 
flavour contribution the sum can be written as F i VFNS (n/ + 1, Q 2 ), provided 
Q 2 ^> m 2 (see Eq. (2.30)). This procedure imposes relations between the 
parton densities appearing in Fi(rif,Q 2 ) and the new ones showing up in 
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Fj VF (rif + 1,Q 2 ). These relations, given by (2.37)-(2.41), have to be satis- 
fied in the VFNS approach. Unfortunately the existing parton-density sets 
in the literature do not satisfy these requirements. They do however satisfy 
the RGE's in (2.52)-(2.55). The consequences of the relations between the 
new parton densities taken at rif + 1 flavours and the old ones presented at rif 
flavours will be discussed for charm production (rif — 3) in the next section. 
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3 Validity of FOPT and the VFNS 



In this section we apply the findings of the last section to charm produc- 
tion. Here we want to investigate which of the two approaches i.e., FOPT 
or VFNS, is the most appropriate to describe the total structure function 
Fi(x,Q 2 ) and its charm component Fi^ c (x,Q 2 ,m 2 ) in the different kinemat- 
ical regimes. To do this we will study the following issues. The first one 
concerns the question at which values of x and Q 2 the large logarithmic 
terms ln l (Q 2 / m c) ln-'(/i 2 /m 2 ) in the heavy-quark coefficient functions consti- 
tute the bulk of the radiative corrections to the charm component of the 
structure functions. Further in the spirit of VFNS one should use the parton 
densities defined in (2.37)-(2.41) instead of the usual ones which are available 
in the literature. In the former there exist direct relations between the parton 
densities taken at rif — 3 (no charm-quark density) and the parton densities 
at n/ = 4 (charm-quark density included). In the literature this relation 
is broken because in all parametrizations the charm-quark density is given 
by an arbitrary function in which the parameters are fitted to the data and 
there exists no direct relation between the charm density and the remaining 
u,d,s and g densities except that they have to satisfy the momentum sum 
rule. The effect of this difference on the charm structure function F 2 ^ c (x, Q 2 ) 
will be investigated. Finally we study the difference in the behaviour of this 
function which occurs when going from the FOPT to the VFNS descriptions. 
In particular we focus on the visibility of the charm threshold which appears 
in FOPT but is absent in the case of the VFNS. We also study the rate of 
convergence of the perturbation series over a wide range of Q 2 - values which 
will be different for these two approaches. 

Before presenting our results we want to mention that all perturbative 
quantities like the operator-matrix elements (OME's), heavy-quark coeffi- 
cient functions and light-parton coefficient functions are presented in the 
MS-scheme. Therefore we have to use parton densities parametrized in the 
same scheme. Furthermore in FOPT the number of light flavours in the 
running coupling constant and the coefficient functions has to be equal to 
three (A 3 = 232 MeV (LO); A 3 = 248 MeV (NLO)) whereas in the VFNS 
this number should be equal to four (A 4 = 200 MeV (LO and NLO)). For 
the mass of the charm quark and the factorization (renormalization) scale 
we have chosen m c = 1.5 GeV/c and jj 2 = Q 2 respectively because the same 
scale is usually adopted for the light-parton components of the structure 
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functions. Notice that in the literature @, |J, [0 a different scale is chosen 
for the charm structure function in FOPT. However the results in NLO are 
rather independent of the scale as is shown in [IT|,[I2] [TH . 

Starting with the first question about the dominance of the large loga- 
rithms we define the charm component in the FOPT approach by the function 
Fi,c{x,Q 2 ,fn 2 ). The latter is computed from (2.29) and contains the heavy- 
quark coefficient functions and H it k, where i — 2, L and k = q,g. The 
above coefficient functions are exactly calculated up to 0(a 2 ) in 0. Their 
asymptotic expressions, which contain the large logarithms above, as well as 



constant terms, are presented in [ID]. The latter are strictly speaking only 
valid when Q 2 3> m 2 . Further let us denote by _p?^ act and F^ mp the charm 
structure functions computed by using the exact and asymptotic forms of 
the heavy-quark coefficient functions respectively. In order to determine the 
Q 2 value above which F^ act and F^/ mp coincide we plot the following ratios 



2 _ F^(x,Q 2 ,m 2 c ) 
i^ act (x,Q 2 ,m 2 ) 

For these plots we adopt the GRV94HO parton density set flT7| . The reason 
that this set is chosen is because it is obtained from a fit to the deep-inelastic 
scattering data performed in the spirit of the FOPT approach, which means 
that the number of active flavours is chosen to be three and the charm com- 
ponent of the structure function is calculated from the photon-gluon fusion 
process and its higher order QCD corrections. 

In Fig. 5 we plot R 2 (NLO) as a function of Q 2 for four different values 
of x i.e. x = 10 _1 , 10 -2 , 10~ 3 and 10~ 4 . From this figure we infer that R 2 
tends to unity at Q 2 values which are an order of magnitude larger than m 2 
( 1 > R 2 > 0.9 for Q 2 > 20 (GeV/c) 2 ). This holds provided x < 10~ 2 , 
where the limit R 2 = 1 is almost always approached from below. However 
for x > lO^ 1 the limit is reached from above and at a much bigger value of 
Q 2 (1.1 > R 2 > 1 for Q 2 > 300 (GeV/c) 2 ), which is two orders of magnitude 
larger than m 2 . The fact that the rate of convergence to R 2 = 1 is much 
slower at large x is more clearly shown in Fig. 6. Here we plot R 2 (NLO) as 
a function of x at Q 2 = 10, 50 and 100 (GeV/c) 2 . At very large x all curves 
strongly deviate from R 2 = 1. When Q 2 decreases the deviations occur at 
smaller x. This fact that R 2 ^ 1 at large x and small Q 2 can be wholly 
attributed to threshold terms, which are present in the exact heavy-quark 
coefficient functions, but which are absent in their asymptotic expressions. 



Mx,Q 2 ,m 2 c )= * 'J - (3.1) 
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In Fig. 7 we show the same plot for Rl(NLO) as presented for R2(NL0) 
in Fig. 5. Here the approach to Rl = 1 starts at a much larger value of Q 2 
than has been observed for R 2 . One sees that 0.9 < Rl < 1.1 when Q 2 > 10 3 
(GeV/c) 2 . Again the rate of convergence is slower as x increases, which is the 
region where threshold effects become important. This becomes even more 
visible in Fig. 8. The reason that threshold effects are dominant at large x 
and small Q 2 can be explained when one looks at the convolution 



dz tx 2 \ TT ( Q 2 m 2 c \ 

L tMp^M* ;?■■?)■ (3 - 2) 



with a similar expression when is replaced by L^. Here and Hi^L^ 
denote the parton densities and the heavy-quark coefficient functions respec- 
tively. The threshold value is given by 

Zth = Q 2 + 4 m 2 c (3 ' 3) 

in the expression for the structure function (2.29). From the above equations 
one infers that when x is very large and Q 2 is very small x — > z t h so that 
only threshold terms can contribute to the integral (3.2). We also would 
like to comment on the phenomenon that the approach to Ri(NLO) = 1 
is much slower for i = L than for i = 2. It originates from the fact that 
the power of the large logarithms appearing in the heavy-quark coefficient 
functions in the case of i = L is one unit smaller than that for i — 2. This 
phenomenon was also observed for heavy-flavour production in the Drell-Yan 
process [0. It appears that the Q 2 value for which the exact and asymptotic 



expressions of the physical quantities coincide is smaller when the powers of 
the large logarithmic terms increase. Notice that the Born contribution to 
the longitudinal coefficient function does not contain logarithms in the limit 
Q 2 3> ml so that it is independent of Q 2 and m 2 c . In this well as 

in some interference terms in the Drell-Yan process, the convergence to the 
asymptotic expressions takes place at an extremely large value of Q 2 . 

We have also studied Ri(LO) in the Born approximation to the charm 
structure functions in (3.1). Here it turns out that Ri(LO) > 1 for all x 
and Q 2 so that the limit is always approached from above. This behaviour 
is different from the one observed in NLO (see Figs. 5-8). Further Ri(LO) is 
closer to unity at small Q 2 and large x than in the case of Ri(NLO) which 
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implies that threshold effects in leading order are smaller than in next-to- 
leading order. However the Q 2 value at which Ri becomes equal to one is 
essentially the same in LO and NLO. In Section 5 of ]TT| we made the same 



study of Ri in (3.1) but on the level of the heavy-quark coefficient functions 
themselves. A comparison with the results from |10j reveals that the Q 2 



value for which the asymptotic and exact heavy-quark coefficient functions 
coincide is the same as the one obtained for the charm structure functions in 
(3.1). Hence the convolution with the parton densities in (3.2) hardly affects 
the Q 2 value at which the exact and asymptotic expressions coincide. 

Next we discuss the parton densities which emerge from the VFNS ac- 
cording to Eqs. (2.37)-(2.41). The most interesting among them is the 
charm-quark density which appears in the four-flavour scheme. It is derived 
from the formula in Eq. (2.38) where we choose nf = 3. Up to 0(a 2 ) Eq. 
(2.38) becomes equal to 

/ c v ™ s (4, x, /i 2 ) ee / 4 (4, x, /x 2 ) + A(4, *, /i 2 ) 
>a s (n 2 )\2 /-i dz 



(3.4) 

where A F ^ q {2 \ A%f(i = 1,2) with H = c are the OME's presented in (B.l)- 
(B.3) which are renormalized in such a way that a s in (3.4) depends on four 
flavours (A 4 = 200 MeV). Further one has to put rif = 3 in the OME's 
above which up to order a 2 are independent of the number of flavours. The 
quantities £(3,/i 2 ) and G(3,fi 2 ) represent the singlet combination of par- 
ton densities and the gluon density in the three-flavour scheme respectively. 
Notice that the 0(a s ) term was already introduced in Eq. (10) of 0. 

In order to make a comparison with the existing charm-quark densities 
present in the literature we now choose the GRV92HO set |19[ to compute 
S(3,yU 2 ) and G(3,/i 2 ) in (3.4). The reason is that the GRV92 set contains 
a charm-quark density which is not included in the GRV94 set |T7||. For a 



comparison between the charm-quark density / c Yc and the one presented by 
the GRV92 set we have to impose the same boundary conditions. Therefore 
we require /^S (4, x, m 2 ) = (3.4) since the same has been done for the 
charm-quark density in fl§|| . To that order the non- logarithmic terms in 
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A^q'^ 1 (B.l) and A^ 2 ^ (B.3) have to be removed, which is not needed for 
A^ 1 ' (B.2) because the latter already vanishes at /i 2 = m 2 . 
In Fig. 9 we plot the ratio 



fc+c (4, 3?, 



-Rch — ,pnp/, ->\ ) (3-5) 



as a function of fi 2 for the values x = 10 _1 , 10 -2 , 10~ 3 and 10" 4 . Here f^™ s 
is computed from (3.4) by choosing one of the parton density sets in the liter- 
ature (here GRV92HO) for the determination of £(3,/i 2 ) and G(3,/z 2 ) while 
ffPjF is the charm-quark density belonging to the same set. For x < 10~ 2 it 
turns out that 1 > i? c h > 0.9. This is surprising because (3.4) is calculated 
up to finite order in perturbation theory whereas the GRV density is the 
NLO solution of RGE's in which all leading and next-to-leading logarithms 
are resummed. Hence we should expect a larger difference between the scale 
evolutions. However from the results in [pOR , pl| we anticipate that this dif- 



ference is not so dramatic as long as /i 2 /m 2 < 10 3 . Apparently the leading 
logarithms beyond 0(a 2 ), which are neglected in (3.4), but not in the GRV 
charm density, do not play an important role provided n 2 is not chosen to 
be too large. What is really striking is that the charm density depends on 
four-light-parton densities in the VFNS, whereas in the case of GRV this 
density is unconstrained except that it has to satisfy the momentum sum 
rule. Therefore it is unexpected that the ratio i? c h is so close to unity. 

Besides the appearance of the charm-quark density, the four other light- 
parton densities representing u,d,s and g are modified while going from three 
to four active flavours. The modification for the light-quark densities is given 
by Eq. (2.37). Summing the latter over the three light flavours (u,d and s) 
one obtains the singlet density corrected up to 0(a 2 ) which is given by 

1 dz. 



X'(W 2 )= / -E(3,-,/x 

Jx Z v Z 



where A^'^ is given in (B.4) with H = c. Notice that S(4, x,/i 2 ) = 
E'(4,x,fi 2 ) + f c+5 (4:,x,fi 2 ) in Eq. (2.39). Up to 0(a 2 s ) the gluon density 
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in the four-flavour scheme (see Eq. (2.41)) is 



,(2)(, V_ 



^S,(2)/ ifV 

4tt ' "' c V ' m 2A 



1 c22 

z 



+ /"s(/i 2 



(3.7) 



S (2) 

The functions and (i = 1, 2) with = c are presented in Eqs. 



|S,(i) 



B.5)-(B.7). Notice that up to 0(a 2 ) the above OME's are independent of 
the number of internal flavours rif. 

In the literature all parton density sets are presented either in the three- 
flavour scheme (n f = 3) e.g. GRV94 JTJ, GRV92 |05| or in the four-flavour 
scheme (n/ = 4) e.g. CTEQ |22] and MRS [23]. This implies that nf, 
which occurs in the anomalous dimensions and in the coefficient functions, 
is kept fixed above and below the charm threshold. In the case of GRV92 
it means that even when Q 2 ^> m 2 the light-parton densities evolve accord- 
ing to rif = 3 whereas in principle nf = 4 has to be chosen. In Fig. 10 
we show the singlet combinations of parton densities E(3,/x 2 ) and £'(4,/i 2 ) 
(3.6) and, as the difference between them is 0(a 2 ), it is essentially invisible. 
Therefore the error made in choosing £(3,// 2 ) instead of E'(4,/x 2 ) above the 
charm threshold is extremely small. In Fig. 11 we make the same study for 
the gluon density where we compare G(3,/i 2 ) with G(4, /z 2 ). Here the latter 
function yields the higher curves. In this case the difference is a larger than 
for the singlet-quark combination (3.6). This is no surprise because the dif- 
ference already starts in 0(a s ) in the gluon case (see (3.7)). We also made a 
comparison with the parton densities of [^2| (CTEQ) and [^3| (MRS). In this 
case the parton-density sets are presented in the four-flavour scheme. Invert- 
ing Eqs. (3.6) and (3.7) one can easily get the corresponding densities for the 
three-flavour scheme and recompute f^™ s in (3.4). Like in the GRV92 case 
the fitted charm-quark density f^+z and the computed charm-quark density 
f^™ s show a remarkable agreement for each parton-density set separately. 
Since f^™ s is mainly determined by the gluon density in (3.4) the mutual 
differences between the charm-quark densities in the various parton-density 
sets may be attributed to the different parametrizations for the gluon. 
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Summarizing our findings above we can conclude that for the actual com- 
putation of the charm structure function F 2 ™ S in (2.60), it does not makes 
much difference when the light-flavour densities in the the VFNS, as given 
by Eqs. (3. 4), (3. 6) and (3.7) are replaced by those obtained from the various 
parton-density sets available in the literature. 

Now we want to study the differences between the FOPT and VFNS 
approaches for the description of the charm structure function F 2jC . Using 
the set GRV92HO @ we present plots for F 2 e * act (2.29) in Figs. 12a-e (see 
also above (3.1)), which represents the FOPT. We do the same for F 2 VFNS , 
which is given by Eq. (2.60), using the densities in Eqs. (3. 4), (3. 6) and 
(3.7). Both are calculated in NLO. For Q 2 = 3 (GeV/c) 2 (see Fig. 12a) we 
observe the threshold effect in _F|* act at large x. Here this function vanishes 
for x > 0.1. This effect does not show up in F 2 VFNS although the latter IS so 
small in the large x- region that, in view of the low statistics of the data in 
that region, any distinction between the FOPT and the VFNS approaches 
will be invisible. The difference becomes even smaller when Q 2 increases 
(see Figs. 12b-e). However at small x the difference between _F|* act and 
F 2 Y becomes more conspicuous. It is in this region that data will appear 
from the HERA collider experiments. Although the difference is not that 
large at Q 2 = 3 (GeV/c) 2 it becomes big for 5 < Q 2 < 10 (GeV/c) 2 (see 
Figs. 12b, c). Here F 2 VFNS exceeds Ff^ ct by more than 60% with respect 
to the latter quantity at x — 10 -4 . The difference becomes less when Q 2 
increases, e.g. at Q 2 = 100 (GeV/c) 2 (Fig. 12e) it is about 25%. The Hl- 
collaboration has published data in [H] for Q 2 = 12, 25 and 45 (GeV/c) 2 which 
lie in the region x < 10 -2 . From their analysis we infer that the predictions 
from FOPT lie below the data which means that, in view of our findings 
above, the predictions from VFNS are in better agreement with experiment. 
However one has to be cautious to draw too premature conclusions from 
the observations made above. First one has to bear in mind that F 2 VFNS is 
derived from F^J mp via mass factorization (see Section 2). Therefore the 
former is only reliable when the large logarithms mentioned in the beginning 
of this section dominate the corrections to the charm structure function. 
From Figs. 5,6 we infer that this happens for Q 2 > 20 (GeV/c) 2 . Another 
criterion is the perturbative stability of the various approaches, which will 
be discussed shortly, and from which we infer that the VFNS is better than 
FOPT for Q 2 > 10 (GeV/c) 2 . This implies that F 2 V C FNS gives a more reliable 
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description of the data than F 2 e * act for Q 2 > 20 (GeV/c) 2 . 

Further we have also plotted F ai ^ ymp (see above (3.1)). As expected from 
Figs. 5,6 this function approaches F|^ act when Q 2 gets very large. We did 
not plot F 2 a ^ ymp at Q 2 = 3 (GeV/c) 2 (Fig. 12a) because here it becomes 
negative for x < 0.1. From Figs. 12b-e one observes a noticeable difference 
between F%*f mi? and F 2 VFNS , which is remarkable since the latter is derived 
from the former via mass factorization; a procedure which holds in all orders 
of perturbation theory. In our actual computations F^J mp has been com- 
puted up to order a 2 whereas F 2 VFNS involves contributions. The latter 
are due to the convolutions in Eq. (2.60) of the order a 2 corrected densities 
in (3. 4), (3. 6) and (3.7) with the order a s corrected light-parton coefficient 
functions. Due to the large logarithms these order contributions are quite 
appreciable and they survive in the region 3 < Q 2 < 100 (GeV/c) 2 covered 
by our plots. The corrections beyond order a 2 present in F 2 VFNS can be re- 
summed using the RGE's in (2.52)-(2.55). This procedure has been applied 
to obtain the parton densities in the literature. If we substitute one of them 
(here we use GRV92HO |19|) in the expression for F 2 VFNS in Eq. (2.60) we 



call the corresponding charm structure function -F 2 P ,P F . The results for this 
function are given in Figs. 12a-e. As can be expected from our discussions 
above F^ F is slightly larger than . This is not surprising because 

the parton densities represent the resummation of the large logarithms in all 
orders of perturbation theory. 

We also want to comment on the presentation of the charm structure 
function in (2.66) originating from (9) in 0] (ACOT). From (2.60) and the 
definitions for F i e £ act ,i^ ymp given above (3.1) one can easily derive that Eq. 
(2.66) is nothing but the 0(a s ) approximation to the following expression 



F t f°^x,Q 2 ,m 2 c ) = F^(x,Q 2 ) + F^(x,Q 2 ,ml) 

-Ft7 mp (^Q 2 ,m 2 c ). (3.8) 



Up to 0(a s ), F^™ s and F^ act — F^ ymp correspond with the first and the last 
part of (2.66) respectively (for the last part see also (2.67)). In the previous 
figures we have seen that for Q 2 > 20(GeV/c) 2 , i^* act approaches F^ ymp so 
that F/^ COT coincides with F^, FNS . Therefore Ff^ OT gives a good descrip- 
tion of the charm structure function at large Q 2 -values. However at small 

Q2 ? F VFNS ^symp ^ p igg gQ ^ ^ACOT ^ nQt dominated by 

^exact Therefore Ff^ OT (3.8) does not have the correct threshold behaviour 
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exhibited by F^ act . The inequality between F^™ s and F^ mp is due to 
the large value of a s {Q 2 ) and the different ways that the corrections beyond 
0(a s ) have been included in the latter two functions. In our case i^ ymp is 
computed up 0(a 2 ) whereas F^™ s contains even O(a^) contributions. If 
F^™ s is replaced by F^ F the above inequality becomes even larger since 
the latter includes all leading contributions in all orders of a s . Due to these 
considerations and the fact that the large logarithms only show up at large 
Q 2 the expressions for -F^, FNS and F^/ mp have no physical meaning in the 
threshold region (low Q 2 ). Hence they should be dropped in the latter region 
so that here the charm structure function is only represented by F^ act . In 
fact the EMC data on charm production were recently reexamined by [p4 



using the NLO results for F|* act (x, Q 2 , m 2 ). The theoretical results are in 
excellent agreement with the experimental data, except for one data point. 

We have also studied the total structure function F 2 (x,Q 2 ), where the 
charm component is included. Since this structure function is dominated 
by the light-parton (u,d,s and g) contributions the differences between the 
various descriptions is much smaller than those observed for the charm com- 
ponent. At maximum these differences are of the order of 10%, which occurs 
in the region 5 < Q 2 < 10 (GeV/c) 2 . 

To study the stability of the perturbation series for F^ c one can proceed in 
two different ways. The first one is discussed in |1 and concerns the behaviour 
of the charm structure function with respect to variations in the factorization 
scale. It was found that near threshold (large x and small Q 2 ) Fff ct shows 
a better stability than F 2 ™ s under variations of the factorization scale. 
Far above threshold it turns out that just the opposite happens, so that 
at large Q 2 it is more preferable to use -F 2 V C FNS instead of F 2 e * act (FOPT). 
However in the analysis of || the NLO corrections from || were not taken 
into account. In ||rTJ,|[12[ and |13| these corrections were included and one 



could show that far away from threshold _p™ act is as stable as F 2 VFNS with 
respect to variations in the factorization scale. The second way to study the 
stability of the perturbation series is to look at the actual size of the higher 
order corrections. They have to decrease when the order in a s increases. To 
be more specific we study the following quantities 



R (i) _ F t ( ! ] c {x,Q 2 ) 
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where denotes the 0(a l s ) corrected charm structure function. Our crite- 
rion is that the perturbation series gets more stable if Kf] -> 1 for increasing 
/. Here we want to compare K^l ,eKact (FOPT) with ^' VFNS (/ = 1,2), which 
are derived from _F|* act and F^™ s respectively. In virtue of the observations 
made in Figs. 12a-e we replace the latter by i r 2 p ( P F because this does not 
appreciably change the results which we now present. In Figs. 13a-g we plot 
i^i' 6 ^ *! -^2!c' PDF an d -^2^c' PDF - Here ^W' exact j s given by the ratio of the 
NLO over the LO charm structure functions in FOPT using the GRV92-set 
[p7gf| . The same holds for K^l' PDF which is obtained from the NLO over LO 

approximations to Eq. (2.60). Finally K^l' FBF stands for the ratio of the 
next-to-next-to-leading-order (NNLO) and NLO corrected charm structure 



functions, where in NNLO only the light-parton coefficient functions of |L6 



have been included. Notice that the NNLO parton densities are not known 
yet because the three-loop splitting functions (anomalous dimensions) have 
not been calculated. From Fig. 13a we infer that for Q 2 = 3 (GeV/c) 2 , 
i^2^' exact is rather close to unity over the whole x-range contrary to K^ff^ 
(I = 1,2) which deviate from unity in spectacular ways. This shows that 
one should use FOPT instead of VFNS at small Q 2 . At Q 2 = 5 (GeV/c) 2 
(see Fig. 13b) the deviation from unity becomes the same for #«> exact a nd 

ifj 1 ] ' PDF in the small ^-region. However at large x (x > 0.01), which repre- 
sents the threshold region, the description by FOPT is still superior to the 
one given by VFNS. This picture changes when Q 2 > 10 (GeV/c) 2 (Fig. 13c). 
In the threshold region {x > 0.01) both approaches are equally bad whereas 
at small x the VFNS is better than FOPT. The superiority of the VFNS 
with respect to the FOPT approach becomes even clearer when one looks 
at Eq. (3.9) plotted at still larger Q 2 -values, which are probably not realiz- 
able experimentally, in Figs. 13d-g. This holds for all x-values including the 
threshold region. We also notice a considerable improvement when the order 
a 2 corrections are included in .F 2 VFNS . In particular K^l^^ is much closer 

to one than K^l^^ at larger Q 2 -values, indicating the rapid convergence of 
the perturbation series in case of VFNS. From the observations made above 
we conclude that Eq. (3.9) provides us with a better criterion to determine 
the predictive power of the perturbation series than the investigation of the 
stability of the charm structure function under changes in the factorization 
scale. It shows very clearly that for Q 2 > 10 (GeV/c) 2 it is better to use 
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F 2 VFNS instead of _F|* act even when x gets large. 

The most important results found in this paper can be summarized as 
follows. When the charm structure function F 2)C is computed in FOPT for 
Q 2 > 20 (GeV/c) 2 the results obtained from the exact (F 2 e * act ) and the 
asymptotic heavy-quark coefficient functions (F 2 a ^, ymp ) are indistinguishable. 
For Fl jC the minimal value of Q 2 becomes larger, say around Q 2 = 10 3 
(GeV/c) 2 . Above this minimal value the large logarithms dominate the per- 
turbation series and they can be resummed after having performed mass fac- 
torization on the heavy-quark coefficient functions. In this way starting from 
^asymp ^ three-flavour scheme one can derive an expression for F 2 V C FNS 
valid in the four-flavour scheme. This procedure imposes a relation between 
the parton density sets parametrized at three and four flavours, which has 
to be satisfied in the VFNS approach. In the literature the parton density 
sets (PDF's) do not obey this requirement although in practice this has no 
serious consequences for the prediction of F 2 VFNS , which is almost equal to 
F 2 P C DF . Further it turns out that F 2 V C FNS > F 2 e * act for all Q 2 - values. This can 
be attributed to the higher-order corrections appearing beyond FOPT which 
are included in VFNS. Finally we have shown that comparisons between the 
charm structure functions calculated in different orders in a s give better in- 
dications about the stability of the perturbation series than a variation in the 
mass-factorization scale. It turns out that below Q 2 = 10 (GeV/c) 2 , _F|* act 
(FOPT) gives a better description of the charm structure function whereas 
above this value it is much better to use -F 2 VFNS . In particular this holds in 
the small x-region under investigation by the HERA experiments. 

We stress again that the VFNS is only valid for totally inclusive quan- 
tities. At the level of exclusive distributions the large logarithms we have 
been discussing simply do not exist. If the experimental analysis is carried 
out on the basis of the photon-gluon fusion model in NLO with three light 
flavours then we expect that the observed charm contribution F 2;C (x, Q 2 , m 2 ) 
should agree with the NLO results F 2 sl ^ ct (x, Q 2 , m 2 ) for small Q 2 (say Q 2 = 5 
(GeV/c) 2 ) and with F 2 v c FNS (x, Q 2 , m 2 ) for large Q 2 (say Q 2 = 100 (GeV/c) 2 ). 
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Appendix A 



In this appendix we present the unrenormalized operator-matrix elements 
^ g q%, A g gjj and h where denote the coefficients of (a s /47r) 1 in 

the perturbation series of the operator-matrix elements (OME's). The cor- 
responding two-loop Feynman graphs are presented in Fig. 3 and Fig. 4 
respectively and the calculation proceeds in the same way as is outlined for 



the other OME's in ||10|| . Using n-dimensional regularization for the ultravi- 
olet and collinear divergences the results are given by (n = 4 + e) 
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Here S e denotes the spherical factor which is given by 

S e = exp{|(7ij-ln47r)}, 
where je is the Euler constant. 
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The last term in Eq. (A. 4), which is proportional to T 2 , is due to the one-loop 



correction to A^ 1 ^ in Eq. (A. 3). This correction is represented by the heavy 
quark (/) loop contribution to the gluon self energy where / represents all 
heavy flavours starting with the quark H (m# = m) and ending with the top 
quark t. The corresponding graph is not shown in Fig. 4. In the next section 
the renormalization of the OME A s ggH will be chosen in such a way that the 
heavy quarks with rrif > m decouple from the running coupling constant. 
This implies that the contributions in the sum of Eq. (A. 4) with / > H 
completely vanish in the renormalized OME A^ g H presented in the next 
section. However the contribution due to / = H remains in the renormalized 
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expressions like those coming from the rif light flavours. This renormalization 
prescription implies that the running coupling constant is presented in the 
MS-scheme and it depends on rif + 1 light flavours including the heavy quark 
H. 

The 1/(1 — z)+ terms appearing in Eqs. (A. 4) and (A. 5) have to be 
understood as distributions, namely 

f dz(-±-) f(z) = I' dz-*—[f{z) - /(l)] . (A.6) 
Jo \1 — z J Jo 1 — z 

The colour factors in SU(N) are given by 



N 2 - 1 1 

C F = C A = N T f = -, (A.7) 

2N 1 2 



with N = 3 for QCD. 
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Appendix B 



Here we present the renormalized operator-matrix elements (OME's) corre- 
sponding to the unrenormalized expressions given in Appendix C of [llj and 
in Appendix A of this paper. All OME'S have been renormalized in the MS- 
scheme. In particular the renormalized coupling constant is presented in the 
above scheme for rif + 1 light flavours. Here the heavy quark H is treated 
on the same footing as the light flavours and it is not decoupled from the 
running coupling constant in the VFNS approach. The (a s /47r) 2 coefficient 
in the heavy-quark OME A^ q is given by 
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The (a s /47r) and the (a s /47r) 2 coefficients of the heavy quark OME A H 
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respectively. Now we present the renormalized expressions for the heavy- 
quark loop contributions to the light-parton OME's denoted by A^h- The 
coefficients of the {a s /Arc) 2 terms in A qqH and A gq ^ H are 



4 NS,(2) / m 
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i r 
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(B.4) 



and 



A 



s,(2) / m 



<iq,h y 
+ 
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(B.5) 



respectively. The coefficients of the a s /47r and (a s /47r) 2 terms in are 



jS,(l) 
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and 
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2 m" 
fi 2 



+ \C F T f 
80 



16 



8(1 + z) In 2 z + (24 + 40z) lnz + 64- 32^ 

3z 



y ^ 2 + 45(1 - z) 



+ C A T f 
184 



16,., Nl 80/1 



+ 



184 232 152 184 2 16 r/1 

H 1 z z 2 H 5(1 - z) 

9z 9 9 9 3 v ; 



ln 



+C F 7) j -(1 + z) In 3 ^ + (6 + lOz) In 2 z + (32 + 48z) ln z 



43 
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+ 



27 



/ I \ 556 
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(B.7) 



[z] and the Nielsen functions 



respectively. 

The definitions for the polylogarithms Li r 
Sn,p( z )-> which appear in the above expressions, can be found in [p5|] . We 
have checked that the renormalized OME's given above satisfy the sum rules 
presented in Eqs. (2.43) and (2.44). This provides us with a strong check on 
our results in [TO] and in this paper. 
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Figure Captions 



Fig. 1. 0(a 2 ) contributions to the purely-singlet parton structure function 
Tfq representing the subprocess 7* + q — > q + q' + q'. Here q and q' 
are represented by the dashed and solid lines respectively. In the case 
of heavy-quark production q' = H and these graphs contribute to the 
heavy-quark coefficient function Hfg. 

Fig. 2. 0(a 2 ) contributions to the purely-singlet OME A^, . Here q and q' 
are represented by the dashed and solid lines respectively. In the case 
of q' = H these graphs contribute to the heavy-quark OME Ajj q . 

Fig. 3. Two-loop contribution to the OME A^ H . The dashed and solid 
lines represent the light quark q and the heavy quark H respectively. 

Fig. 4. Two-loop graphs contributing to the OME A^ g H . The dashed and 
solid lines represent the Faddeev-Popov ghost and the heavy quark H 
respectively. The graph with the external Faddeev-Popov ghost (fig.4g) 
has to be included if the sum over the gluon polarization states involves 
the contributions from unphysical polarizations. 

Fig. 5. R,2(NLO) (3.1) plotted as a function of Q 2 at fixed x; x = 1CT 1 
(dashed-dotted line), x = 1CT 2 (dotted line), x = 10~ 3 (dashed line) 
and x = 10~ 4 (solid line). 

Fig. 6. R 2 {NLO) (3.1) plotted as a function of x at fixed Q 2 ; Q 2 = 10 
(GeV/c) 2 (dotted line), Q 2 = 50 (GeV/c) 2 (dashed line), Q 2 = 100 
(GeV/c) 2 (solid line). 

Fig. 7. Same as in Fig. 5 but now for Rl(NLO). 

Fig. 8. Same as in Fig. 6 but now for Rl(NLO). 

Fig. 9. The ratio i? c h = /^ NS //™ F (3.5) as a function of \j? at fixed x 
in NLO. The charm densities fj™ s and are given by (3.4) and 

the set GRV92HO |]| respectively; x = 10" 1 (dashed-dotted line), 
x = 10~ 2 (dotted line), x = 10~ 3 (dashed line), x = 10~ 4 (solid line). 
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Fig. 10. E(3, x, fi 2 ) and S'(4, x, /i 2 ) (the difference between both singlet 
combinations of quark densities is unnoticeable in the figure) plotted 
as functions of /i 2 at fixed x; x = 1CT 1 (solid line), x = 1CT 2 (dashed 
line), x = 1CT 3 (dotted line), x = 1CT 4 (dashed-dotted line). 

Fig. 11. G(3,x,/i 2 ) (lower lines) and G(4, x, /i 2 ) (upper lines) plotted as 
functions of fi 2 at fixed x; x = 1CT 1 (dashed-dotted line), x = 1CT 2 
(dotted line), x = 1CT 3 (dashed line), x = 1CT 4 (solid line). 

Fig. 12a. The charm component of the structure function given by Fi, c in 
NLO as a function of x at Q 2 = 3 (GeV/c) 2 ; F 2 cxact (x, Q 2 , m 2 ) (solid 
line), F 2 v c FNS (4,x,Q 2 ) (dashed line), F 2 P C DF (4, x, Q 2 ) (dotted line) with 
PDF=GRV92HO. 

Fig. 12b. The charm component of the structure function given by F 2iC in 
NLO as a function of x at Q 2 = 5 (GeV/c) 2 ; F 2 exact (x, Q 2 , m 2 ) (solid 
line), F 2 v c FNS (4,x,Q 2 ) (dashed line), F P C DF (4, x, Q 2 ) (dotted line) with 
PDF=GRV92HO , and F| s c ymp (x, Q 2 , m 2 ) (dashed-dotted line). 

Fig. 12c. Same as in Fig. 12b but now for Q 2 = 10 (GeV/c) 2 . 

Fig. 12d. Same as in Fig. 12b but now for Q 2 = 50 (GeV/c) 2 . 

Fig. 12e. Same as in Fig. 12b but now for Q 2 = 100 (GeV/c) 2 . 

Fig. 13a. The ratios K^\{x, Q 2 ) (eq.(3.9)) plotted as functions of x at Q 2 = 
3 (GeV/c) 2 ; K^'*™* (solid line), K^/ DF (dashed line), K^/ DF (dot- 
ted line), with PDF=GRV92HO. 

Fig. 13b. Same as in Fig. 13a but now for Q 2 = 5 (GeV/c) 2 . 

Fig. 13c. Same as in Fig. 13a but now for Q 2 = 10 (GeV/c) 2 . 

Fig. 13d. Same as in Fig. 13a but now for Q 2 = 50 (GeV/c) 2 . 

Fig. 13e. Same as in Fig. 13a but now for Q 2 = 100 (GeV/c) 2 . 

Fig. 13f. Same as in Fig. 13a but now for Q 2 = 10 3 (GeV/c) 2 . 

Fig. 13g. Same as in Fig. 13a but now for Q 2 = 10 4 (GeV/c) 2 . 
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Fig. 4 ( continued ) 
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